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ABSTRACT 

We present generalized Schiff's transformation on electric dipole moments 
(EDM) in quantum field theory. By the unitary transformation, the time and 
parity violating interaction i^^o ^^F^ is transformed into a new form, but its 
nonrelativistic reduction has a unique form, which is free from Schiff's theorem. 
The relativistic corrections to the new EDM operator turn out to be a small 
increase to the EDM as given by b2(aZ) 2 with f>2 — 2. Therefore, the calculation 
of the EDM with nonrelativistic Hartree-Fock wave functions presents the most 
conservative but reliable estimation for the enhancement factor of the EDM in 
atoms. 
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1. Introduction 

Symmetries in physics are fundamental for understanding nature. The high- 
est symmetry in quantum field theory is the CPT ( charge conjugation, 
parity and time reversal) invariance. Any field theoretical models should 
be consistent with the CPT theorem. The next highest symmetry may be 
the time reversal invariance. Indeed, the T-invariance is also kept well in 
most of the field theoretical models. 

It is of fundamental interest to find the T-violating interactions in nature. 
Electric dipole moments (EDM) of particles, nuclei and atoms in ground 
state reveal the violation of the T-invariance. Until now, the upper limit of 
the neutron EDM is around —(3 ± 5) x 1CT 26 e • cm [1]. 

There have been many experimental efforts to measure the EDM of the 
atomic systems. The best example is found for the EDM of 199 Hg [2]. In 
this case, however, one has to be careful for extracting the EDM of the 
electrons or nucleons from the measurement of the atomic EDM since there 
is Schiff's theorem [3]. This theorem states that the EDM of the atom is 
canceled out due to the symmetry restoration mechanism as long as the 
constituents are interacting through the electromagnetic interactions with 
the nonrelativistic kinematics. 

Due to the presence of Schiff's theorem, people have long believed that the 
calculation of the enhancement factor in atomic systems should be carried 
out by employing the relativistic many body wave functions. However, it is 
not clear to what extent the Dirac Hartree-Fock method for the relativistic 
many body wave functions can be reliable. In fact, the Dirac Hartree-Fock 
method may well have some conceptual difficulties since the relativistic 
many body theory should be treated by field theories. The relativistic 
quantum mechanics is well defined only for the one body case. 
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In this sense, it is of particular importance to find a way to estimate the 
EDM enhancement factor which does not depend on the ambiguity arising 
from atomic wave functions. 

In this paper, we present a generalization of Schiff's transformation which 
is extended to the quantum field theory of electrodynamics. We show that 
the T and P-violating interaction igeipa^^ipF^ is transformed into a new 
form which has a unique shape of the nonrelativistic reduction. For this 
operator, we can reliably calculate the EDM of the atomic system with 
nonrelativistic Hartree-Fock wave functions. 

In particular, the estimation of the EDM with Dirac wave functions in 
hydrogen-like atoms shows that we can write the EDM dz of the atoms 
with Z in units of electric EDM d e as 

d z = ^i(Za) 2 (l + b 2 {Zaf) d e (1.1) 

where b ,bi and e are dimensionless constants which are determined by 
the nonrelativistic wave functions. 62 is a positive number around 62 — 2. 
Eq.(l.l) indicates that the relativistic corrections are of the order of (aZ) 2 , 
and therefore they cannot be very large. Further, the relativistic effects 
must be smaller for the Hartree-Fock case than the hydrogen-like atoms 
since the Hartree-Fock effects always push the electron of the last orbit 
outward and thus the relativity becomes less. 

On the other hand, many body effects on the EDM are very large. This is 
mainly due to the fact that the EDM operator is, roughly speaking, pro- 
portional to ^2, and thus the many body wave functions which are pushed 
outward due to electron electron repulsions reduce the EDM expectation 
values by an order of magnitude compared to the hydrogen-like wave func- 
tion. Thus, the careful treatment of the many body effects is much more 
important than the relativistic corrections to the EDM. 
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In this paper, we calculate the EDM for the hydrogen, Li and Cs atoms 
as examples and estimate the EDM with the nonrelativistic Hartree-Fock 
wave functions. In particular, it is found that the enhancement factor of 
Cs atom in comparison with the electron EDM d e is 91.2, which should be 
compared with the previous estimations of 100-150 [4-10]. 

This paper is organized in the following way. In the next section, we explain 
briefly Schiff's theorem in the nonrelativistic quantum mechanics. Then, 
section 3 treats the field theoretical version of the Schiff transformation. 
In section 4, we evaluate the Foldy-Wouthuysen transformation of the rela- 
tivistic EDM operator and obtain a unique form of the EDM operator which 
has no influence from Schiff's theorem. In section 5, we carry out the nu- 
merical evaluations of the EDM with the hydrogen-like wave functions. In 
particular, we compare the calculations with nonrelativistic and relativistic 
wave functions. In section 6, we present our calculations of the EDM with 
nonrelativistic Hartree-Fock wave functions. Section 7 summarizes what we 
have clarified in this paper. 



2. Schiff's theorem in quantum mechanics 

Here, we first review the shielding mechanism which was initially found by 
Schiff [3] . We consider the hamiltonian Hq r of the atomic system which is 
interacting through the Coulomb force as well as the external electric field 

H o R = E (|f: " eAto) + er t ■ E e A (2.1) 
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where A)( r i) is given as 

MTi) = T-l^wh-\- (2 - 2) 

1 1 * \*-l Lj\ 

To this system, we consider the T-and P-violating interaction H^ d R 

H edm = - ^ • E< (2.3) 

i 

where Ej is defined as 

E i = -V i A (ri)+E ext (2.4) 

dj represents the EDM of the i-th particle and can be written as 

di = gecri (2.5) 

where cr^ represents the spin operator of the i-th particle, g is a constant 
and represents the strength of the T and P-violating interaction, and it has 
the dimension of the length. 

Schiff 's theorem states that the total hamiltonian H NR = Hq R + H^ d R is 
related to the Hq r by the unitary transformation ( Schiff transformation ) 
in the following way 

H NR = H NR + H NR = exp (_ • ^ . ^ R NR ^ ^ £ ^ . ^ 

(2.6a) 

since we can easily evaluate the righthand side of eq.(2.6a) 

exp I p . . d ^ ff^exp (^E^ • d *) = H o R ~ E d * • E* + OQ? 2 ). 

(2.66) 

Therefore, up to the order of g, eq.(2.6a) holds. By the unitary transfor- 
mation, one obtains the same energy spectrum between H NR and Hq R . 
Therefore, this means that the effect of the H^ d R is absorbed into the origi- 
nal hamiltonian and thus cannot be observed at all. This is Schiff 's theorem. 



In his paper, Schiff argued that the shielding of the EDM may be violated 
by the relativistic effects or some other effects like the addition of the strong 
interactions. 

In what follows, we consider the Schiff transformation for the relativistic 
case. 



3. Schiff transformation in QED 

(a) Schiff transformation in lagrangian formulation 

Now, we want to generalize Schiff 's theorem to the field theory. Here, 
we only consider the QED. The lagrangian density for fermions inter- 
acting with the gauge field can be written as 

C = i>{i lv £>» - m )4> - -F^ (3.1) 

where and F^ u denote the covariant derivative and the field tensor, 
respectively and are given as 

D„ = d M + ieA„ (3.2a) 

= d^A v - d v A,. (3.2b) 

In addition to the above lagrangian, we consider the T-and P-violating 
interaction which may presumably be induced by the supersymmetric 
model with some soft breaking interactions [11]. The T-and P-violating 
effective interaction can be given as 

Cedm = -i^^o^vl^F^ (3.3) 
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where a^ u is defined as 

i 

<V= 2^7^-7^)- (3-4) 

The corresponding hamiltonian of eq.(3.3) which depends on the elec- 
tric field can be written as 

HfL = -gejoV ■ E (3.5a) 

where we define the relativistic spin operator £ as £ = 7°7 5 7. 

Here, if we take a simple nonrelativistic limit of this operator, the 
hamiltonian becomes 

H?L ^ -ge<r • E. (3.56) 

Since the d is written as d = geer, this just corresponds to the H^ n . 
Now, we consider the following unitary transformation 

ip' = exp (igi&fSf) ip. (3.6) 

Under this transformation, the total lagrangian density C = C + C e d m 
becomes up to the order of g 

C = C - 2ig4> lb {p^ - eA^. (3.7) 

This shows that the unitary transformation of eq.(3.6) completely can- 
cels out the Cedm term of eq.(3.3), but a new term 2igip^ 5 p^p' J ' / ip ap- 
pears. 

Now, we want to make the corresponding hamiltonian. In this case, 
we should make the conjugate momentum n^ for the A and 11^ for 
the ip. It is clear that the terms proportional to the p disappear from 
the hamiltonian. Therefore, the new EDM hamiltonian becomes 

H^ m = 2 ig ^ 5 (pl-P 2 )- (3-8) 
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This is a new EDM hamiltonian which should be compared with the 
original one in eq.(3.5a). We should also note that the similar EDM 
hamiltonian is obtained in ref.[12]. In this case, the term is missing 
since they treat it in terms of the hamiltonian. Instead, we treat the 
problem in a covariant fashion. However, it is easy to see that the p\ 
term does not contribute to the mixing of the wave functions since it 
does not depend on the coordinates. In this respect, we reproduce the 
result obtained in ref. [12] with the lagrangian formulation with the full 
relativistic covariance. 

The important point is that the new EDM hamiltonian is free from 
Schiff's theorem. Therefore, this term indeed contributes to the gen- 
eration of the EDM. 

Here, we should comment on the physical meaning of the unitary trans- 
formation of eq.(3.6). Under the transformation, we obtain a new la- 
grangian density which does not couple to the external electromagnetic 
field. This is just the statement which corresponds to Schiff's theorem. 
Namely, the EDM can only be induced from the second order effect of 
the new hamiltonian as will be treated later in eq.(5.2). 

(b) Schiff transformation in hamiltonian formulation 

Next, we show here the same type of proof for the Dirac hamiltonian 
since the EDM leftover can be also obtained from the Dirac hamilto- 
nian. This method is similar to the above proof of the field theory, 
but the transformation operator is different from the field theory treat- 
ment, but is the same as the nonrelativistic Schiff transformation. 
The Dirac hamiltonian for fermions interacting with the Coulomb force 
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can be written as 

N 

H£ = £ (Pi • "i + ^7° - eA ( ri ) + er, • E e *') . (3.9) 

i 

In this case, the EDM hamiltonian which corresponds to C e dm can be 
written as 

N 

H™ m = -geJ2lhi-Vi. (3.10) 

i 

Here, we can rewrite the Hamiltonian H^ m as 

N N 

H*L = -ge^i ■ E * - ^E(7- - 1)S 4 • E<. (3.11) 

j j 

In this case, one can easily check that the first term of eq.(3.11) is 
canceled out by the Schiff transformation in the same way as the non- 
relativistic case of eq.(2.6a), and the following equation holds up to 
the order of <?[13]. 

<-^£(7°-l)£i-E, = exp (i^gpi • (H«+H%Je W (-i^gpt ■ 

(3.12) 

Therefore, the leftover is just the second term of eq.(3.11), and we 



define H^ m as 



^ 3 m = -^E( 7 ,°-l)S i -E,. (3.13) 



In the next section, we will prove that the two terms (eq.(3.8) and 
eq.(3.13)) have the identical shape of the nonrelativistic reduction un- 
der the Foldy-Wouthuysen transformation. This nonrelativistic EDM 
operator is free from Schiff 's theorem. 
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4. Nonrelativistic reduction of EDM operator 

When we want to evaluate operators which are expressed in the relativistic 
form, we should prepare wave functions which are obtained relativistically. 
In atomic systems with high Z, the lowest state is indeed quite relativistic. 

However, the system we are treating is not the hydrogen-like atom but 
real atoms. As an example, let us consider the Cs atom. In this case, 
there are many electrons around, and the state for the last electron in the 
ground state is 6s state which is far from relativistic. Therefore, we should 
rather obtain the EDM operator which is reduced to the nonrelativistic 
form. Since the EDM operators obtained in the last section are free from 
Schiff's theorem, we should make the reliable nonrelativistic reduction of 
the operator. 

Here, we employ the Foldy-Wouthuysen transformation which is a unitary 
transformation. It is here noted that the nonrelativistic reduction should 
be done at the last stage. Otherwise, one often makes mistakes since the 
Foldy-Wouthuysen transformation does not necessarily commute with other 
unitary transformation. Since the nonrelativistic reduction is an approxi- 
mation, one has to do it after one has made other unitary transformations. 

In the case of the EDM operator, we have obtained the two different hamil- 
tonians, eqs.(3.8) and (3.13). Here, we rewrite them again as the one body 
operators, 

H™ m = 2ig l0l M-V 2 ) (3.8) 

Hedm — ~9 e (l° — 1)E • E. (3.13) 

They look very much different from each other. Now, we want to make the 
nonrelativistic reduction of the two operators. 

By the Foldy-Wouthuysen transformation which is the most reliable method 
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for the nonrelativistic reduction, we write the reduced hamiltonian after 
some repeated operations of the transformation [14], 

^ = 7.(-+^-^)+£-^[0,[0,£]] (4.1) 

where O and 8 denote the odd and even operators in the gamma matrix 
space. 

In the case of the EDM hamiltonian of eq.(3.8), the O and 8 can be written 

as 

O = a • p - 2ic/7o7 5 p 2 (4.2a) 

S = -eA + er • E ext . (4.26) 

On the other hand, for the EDM hamiltonian of eq.(3.13), the O and 8 
become 

O = a • p (4.3a) 

8 = - e A - #e(7° — 1)S • E + er • E ext . (4.36) 

After some straightforward calculation, we obtain the nonrelativistic hamil- 
tonian for both of the cases of eq.(3.8) and eq.(3.13), 

a ■ E)V 2 - p(r)(V • a) - 2(E • V)(<r • V)l (4.4) 



V - ge 

v edm 



2m 2 

where p(r) and E are defined as 

p(r) = -VMo(r), (4.5a) 

E = -VA)(r) + E exi , (4.56) 

Mr) = — +[ /^Vr' (4.5c) 
r J \r — r'| 

with p (r) = -eE„ \^n{r)\ 2 ■ 

This is a proof that the apparent two different shapes of the relativistic 
hamiltonian eqs.(3.8) and (3.13) reduce to the identical and thus unique 
nonrelativistic hamiltonian. 
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As we stressed before, this hamiltonian is free from Schiff 's theorem. There- 
fore, there is definitely some mixture due to this operator between the op- 
posite parity states like the 6s and 6p states. 



5. EDM with hydrogen-like wave functions 

Since we obtain the EDM operator which is free from Schiff's theorem, 
we can now reliably calculate the EDM in actual cases. Before going to 
the calculations with the Hartree-Fock wave functions (in section 6), we 
first present our calculations with the hydrogen-like wave functions. The 
hydrogen-like wave functions are obviously far from reality in atoms since 
the interactions between electrons are very important, and thus the wave 
functions in the hydrogen-like atoms are quite different from the Hartree- 
Fock wave functions. But the use of the hydrogen-like wave functions helps 
us understand the basic structure of the EDM operator. This is quite 
important since we can evaluate all the matrix elements analytically. In 
addition, it turns out that the EDM eveluated by the hydrogen-like atom 
wave functions presents the upper limit of the enhancement factor even 
though the overestimation is quite significant. 

Here, we should note that the energy eigenvalues are replaced by the esti- 
mation using Hartree-Fock wave function since otherwise it is in fact mean- 
ingless to use the energy eigenvalues of the hydrogen-like atoms which are 
degenerate between nsi and npi states. 

2 2 

Here, we consider the atomic systems in which one electron is found in the 
outer shell like Li or Cs cases. 
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Now, we prepare the nsi and ran state wave functions which can be written 

2 2 



as 



2 



i2n.(r)ei, (5.1a) 

(5.16) 

where R ns (r) and R np (r) represent the radial part of the atomic wave func- 
tions. £i denotes the spin part. 

In this case, we obtain the effective EDM dz for nsi state in atoms as 

2 

< nsi | J2i z i\ n 'pi >< n'pilVedmlnsi > 



"2eE 



— E , 

ns i J-'n'p i 
7 2 



(5.2) 



Here, we only consider the n' — n case since this gives a dominant contri- 
bution to the EDM. Further, we consider the atomic system when only one 
electron is found in the outer shell. Therefore, the summation over % is just 
one single particle only. Also, we neglect the term which depends on the p 
in eq.(4.4) since it vanishes for the hydrogen-like atoms. 

Now, we first evaluate the angular parts of the matrix elements < nsi\z\npi > 

2 2 

and < npi\V e dm\nsi >, and obtain 



< ^ns\z\i(j„p >= - < R ns \r\R np > 



2m 2 



Za 



3 

R np (r) 



1_ (2 dR ns (r) _ d 2 R ns (r) 



dr 



dr 2 



(5.3a) 

r 2 dr. 
(5.36) 



Now, we want to make the matrix elements dimensionless and therefore we 
write them as 



ii ^ h a ° 
< nsi\z\npi >= o — 

2 2 Zj 

< npi\V edm \nsi >= b l7 ^—(Za) (—\ 
2 2 2m 2 Van/ 



(5.4a) 
(5.46) 



where an denotes the Bohr radius in hydrogen atom and is written as 

h 2 



a = 



me" 



(5.5) 
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Further, we write the energy difference AE = E npi — E nsi as 

AE = e m(Za) 2 . (5.6) 
In this case, we can write the EDM as 

d z = ^h(Za) 2 d e . (5.7) 

In what follows, we estimate the values of the &o and b\ which depend on 
the states. Here, before doing so, we should make a comment concerning 
the first order contribution of the V e d m term which couples directly to the 
external field of E ext . This contribution is a few orders of magnitude smaller 
than that of eq.(5.2). However, in the case of hydrogen ground state, the 
contribution amounts to 17 % of eq.(5.2). Therefore, we take into account 
the first order contribution to the hydrogen case. 

(a) Li (Z=3) : nonrelativistic case 

In this case, the electron is in the 2s i orbit. For the bo and b±, we can 
easily evaluate them and find 

b = — v/3 (5.8a) 

h = (5.8.) 

In this case, we can evaluate the EDM assuming the energy difference 
from the estimation using the Hartree-Fock wave functions 

vnot 

AE = 0.135^-. (5.9a) 

Thus, we obtain 

e = 7.5 x 10" 3 . (5.96) 
Therefore, we finally obtain for the EDM value of Li case 

d Li = 0.0324. (5.10) 
14 



This value should be compared with the one obtained by the Hartree- 
Fock calculation. The EDM expectation value with the hydrogen-like 
wave function overestimates the EDM value by a factor of 5 as we will 
see later. 

(b) Li (Z=3) : relativistic case 

Here, we evaluate the EDM matrix element using Dirac wave functions 
which are solved in a pure Coulomb potential V(r) = — ^f-. In this 
case, the wave function is specified by the quantum number n and k. 
The nsi state (k — —1) can be written as 

*n S1 = f f\; Ki } ) • (5.11a) 
Also, the npi state (« = 1) can be written as 

*-j - ( Sp'% ) ■ < 5 - 116 > 

The analytic expressions of the radial wave functions G n>K , F n>lt , G njK 
and F n>K are given in Appendix. 

In this case, we can first evaluate the angular parts of the EDM matrix 
elements as, 

1 f°° 

< * nsi \z\V npi >= -- / (G n _iG n ,i + F^F^ydr (5.12a) 



2 2 



3 Jo 



roc _ 

<^ns 1 \Vedm\^n Pl >= ~2gZa F^F^dr. (5.126) 

2" 2 JO 

Now, for Li case, the electron is in the n = 2 and k — — 1 state. After 
some calculations, we obtain for 2si — 2pi case, 

2 2 

< *2. i H**, >= -VI (l - ^(QZ) 2 ) (I) (5.13a) 

< * 2<} I^|* 2P , >= -£(Z«) (1)' ^ (l + f (aZf 

(5.136) 
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Thus, we obtain the EDM with the relativistic wave functions as 

d£ = 0.032 (l + 1.75(aZ) 2 ) d e . (5.14) 

Therefore, we can see that the relativistic correction to the EDM ex- 
pectation value is of the order of (aZ) 2 , and thus it is quite small in 
this case. 



(c) Cs (Z=55) : nonrelativistic case 

In this case, the electron is in the 6s i orbit. For the &o and b±, we can 

2 

easily evaluate them and find 

b = -3^35 (5.15a) 




(5.156) 



In this case, we can evaluate the EDM assuming the energy difference 
from the Hartree-Fock estimation 

e = 1.31 x 10~ 5 . (5.16) 

Therefore, we finally obtain for the EDM value 

d Cs = 2.66 x 10 3 d e - (5.17) 

This value should be compared with the one obtained by the Hartree- 
Fock calculation. The EDM expectation value with the hydrogen-like 
wave function overestimates the EDM value by a factor of 30 as we 
will see later. 

(d) Cs (Z=55) : relativistic case 

In this case, the electron is in the n = 6 and k = — 1 states. After 
some calculations, we obtain for 6s± — 6»i case. 

2 2 

< ^ 6s , |z|*6p, >= -3V/35 (l - |j(^) 2 ) (f ) (5.18a) 
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2m 2 



9 




(5.186) 



Thus, we obtain the EDM with the relativistic wave functions as 



Therefore, we can see that the relativistic correction to the EDM ex- 
pectation value is again of the order of (aZ) 2 , and thus it is 30 % 
correction to the EDM. But we should note that the relativistic effect 
must be smaller for the relativistic Hartree-Fock case, and thus the 
relativistic correction to the EDM can be at most 30 %. 
Here, we note that the contribution due to the next order term (aZ) 4 
is less than 10 %, and thus we do not have to worry about the higher 
order terms. 

(e) Hydrogen atom 

Finally, we want to present our calculation for the hydrogen atom. In 
this case, there is no stable 2si state, and therefore there is no en- 

2 

hancement for the EDM. However, the hydrogen atom is best studied 
in many respects, and thus there might be some chance that one can 
observe the EDM for this case with very high precision. 

Obviously, we have the wave functions analytically, and thus we can 
calculate the EDM in a closed form. Since the 2si state is not stable, 



we consider the ground state (lsi) which can be mixed with the 2pi 
state. 

For the b and b±, we can easily evaluate them and find 




(5.19) 



2 



>0 



128V6 
729 



(5.20a) 



&i 



16 



(5.206) 



9^6' 
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Since the energy difference is obtained as 



AE = -ma 2 
8 



(5.21a) 



we obtain 



3 



(5.216) 




Finally we obtain the EDM for the hydrogen 



d H = -5.32 x l(T 5 d e . 



(5.22a) 



which should be compared with Sandars calculation [5], 



d H = -2a 2 d e = -10.7 x 1(T 5 4. 



(5.226) 



The difference between them is due to the fact that Sandars took into 
account all of the npi states while we consider only the first excited 
state of 2pi . 

2 

This is quite small and therefore it is not very clear whether one can 
choose the hydrogen atom for the EDM experiment. 



6. EDM from Hartree-Fock wave functions 

In the previous section, we have presented the EDM calculations with 
hydrogen-like wave functions. Now, we evaluate the EDM with the nonrel- 
ativistic Hartree-Fock wave functions [15]. 

Unfortunately, we cannot obtain the Hartree-Fock wave functions analyt- 
ically, and thus we have to carry out all of the calculations numerically. 
However, since the basic expressions are given in terms of the hydrogen-like 
wave functions, we can follow the same notations. 
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(a) Li (Z=3) : Hartree-Fock wave function 

In this case, the electron is in the 2s 1 orbit. For the 6 and 61, we can 

2 

evaluate them and find 

b = -4.17 (5.23a) 

61 = -0.029. (5.236) 

Now the energy difference between 2si and 2pi states is calculated to 



be 



Therefore, we obtain 



Thus, we finally obtain 



,ma 2 



AE = 0.135-—. (5.24a) 



e = 7.5 x 10~ d . (5.246) 



(Iff = 0.00774- (5.25) 

This value should be compared with the EDM value calculated from 
the hydrogen-like wave function 

d Li = 0.0324- 

The absolute value of the EDM is smaller for the Hartree-Fock calcu- 
lation than the hydrogen-like case by a factor of 5. 
This is quite easy to understand since the wave function for the Hartree- 
Fock wave function should be pushed outward compared with the pure 
Coulomb case due to the electron-electron repulsion. 

(b) Cs (Z=55) : Hartree-Fock wave function 

In this case, the electron is in the 6s 1 orbit. For the 6 and b±, we can 

2 

evaluate them and find 

6 = -116 (5.26a) 
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&! = -6.39 x 10~ 5 . (5.266) 
The energy difference between 6si and 6pi states is calculated to be 

o 

AE = 0.025^. (5.27a) 

Thus, we obtain 

e = 1.31 x 1(T 5 . (5.276) 
Therefore, we finally obtain 

dgf = 91.24- (5.28) 

This is the enhancement factor of the EDM in Cs atom from the 
nonrelativistic Hartree-Fock calculation. This should be compared 
with the EDM calculated by the Dirac Hartree-Fock wave function. 
They predict the EDM in Cs atom around (100 — 150)<i e depending on 
the calculations [4-10]. This is slightly larger than the one calculated 
here. This is indeed reasonable since the relativistic effects in EDM 
must be of the order of {aZ) 2 increase compared to the nonrelativistic 
evaluations. 

In any case, it is now clear that the EDM in atoms is indeed enhanced 
compared with the electron EDM d e by a large enhancement factor 
for the atom with relatively high Z. Therefore, it would be better to 
use the atomic systems to observe the EDM. 



7. Conclusions 
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We have presented the nonrelativistic reduction of the EDM operator which 
is free from Schiff's theorem. This is important since we know definitely 
that the effect of the T-and P-violating interaction can be well observed in 
the atomic systems. Further, the enhancement factor is indeed quite large 
for the Cs atom. There are two reasons of the enhancement of the EDM in 
atoms. The first one is that the EDM operator is practically proportional 
to ^ and thus the EDM expectation value become large for the large Z. As 
the second reason, the EDM becomes enhanced due to a very small energy 
difference between the ground state and the first excited state which has 
an opposite parity. Since the Hartree-Fock wave function is quite reliable 
in atomic systems, the estimated value of the EDM with the nonrelativistic 
reduction is indeed reliable. Here, we show that the relativistic correction to 
the EDM expectation value is of the order of iZa) 2 and therefore is not very 
large. Since the relativistic correction tends to increase the absolute value of 
the EDM, the nonrelativistic evaluation gives a conservative enhancement 
factor for the EDM in atoms. 

Here, we should make comments as to what is the difference between the 
present approach and the Schiff's calculation. The basic point is that the 
nonrelativistic reduction should be made at the last stage. For example, 
the simplest EDM operator is obtained when we make the nonrelativistic 
reduction at the lagrangian level. In this case, the EDM operator becomes 
eq.(3.56) as discussed before. This is the worst case since there is no EDM 
operator left due to Schiff's theorem. 

The next level of the approximation which we have not presented in this pa- 
per is that we make the Foldy- Wouthuysen transformation before the Schiff 
transformation. Namely, we make the Foldy- Wouthuysen transformation 
for eq.(3.5a). In this case, we obtain the nonrelativistic EDM operator 
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which is indeed different from eq.(4.4) after we remove the term described 
in eq.(2.3) by the Schiff transformation of eq.(2.66). 

The procedure we employ in this paper is that we make the Schiff trans- 
formation at the lagrangian level and then make the Foldy-Wouthuysen 
transformation to the EDM operator. In this case, it is found that we ob- 
tain the unique nonrelativistic EDM operator which is free from Schiff's 
theorem. 

Therefore, as long as we employ the EDM operator of eq.(4.4), then we can 
estimate the EDM in atomic systems quite reliably with the Hartree-Fock 
wave functions. 

Finally, we make a comment as to how the present calculation is related 
to the relativistic Hartree-Fock calculations. If we want to treat very high 
Z atoms, then we have to carry out the calculations relativistically. In 
this case, the direct evaluation of eq.(3.8) or (3.13) should be done if one 
obtains reliable wave functions with relativistic Hartree-Fock method [16]. 
However, if one really has to treat many body problems relativistically, then 
one has to carry out the calculations field theoretically. But this must be 
an extremely difficult task. 
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Appendix 

Here, we present the radial parts of the Dirac wave functions in a pure 
Coulomb potential V(r) = —^f-- 



1. 2si and 2pi states 

2 2 



<Vi(s) = 



3/2 



2A 5 (l + WQr(2 7 + 2) t /2 

r(2 7 + l)^ Za(Za + A) ' ' 



l-(-l-^) F(-l,27 + l,x)" 



and 

G 2 ,i(z) = 



_ m 3/2 



2A6(l-W)r(2 7 + 2) _ x 



r(2 7 +l)\ Za(Za + A) 



m 



3/2 



2A5(l + WQr(2 7 + 2) t /2 



r(2 7 +l)\ Za(Za-X) 



l-^)F(-l,2 7 +l,x) 



where 



-m 1 



3/2 



2A 5 (1 - W0r(2 7 + 2) 



r(2 7 + l)\ Za(Za-A) 



7 



= V 1 - (^«) S 



A= Wl- 



E 



rn 



E 1 
W = — = 



m /i+W 



x = 2mAr 



Here, the function F(a, b, c) denotes the confluent hypergeometric function. 
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2. 6si and 6m states 

2 2 



G 6 „i(x) 



3/2 



2A*(l + W0r(2 7 + 6) , , x 



r(2 7 + l)^ 5!Za(Za + A) 



X 



-l--^jF(-5,2 7 +l, 



_ m 3/2 



2A 5 (1 - W^)r(2 7 + 6) 



r(2 7 + l)^ 5!Z«(Za + A) 



and 



x 



5F(-4, 2 7 + 1, x) - (-1 - -^J F(-5, 2 7 + 1, x) 



G 6 ,i(x) 



rn 



3/2 



2A5(l + ^)r(2 7 + 6) , /2 



r(2 7 + l)\ 5!Za(Za-A) 



x 



F(-5,2 7 + l,x) 



^8,1 (*) 



-m 



3/2 



2AB(l-W)r(2 7 + 6) _ x 



r(2 7 + l)\ 5\Za(Za-\) 



where 



x 



(Z cx \ 
1 - F(-5,2 7 + l,x 



7 = Ji-(z a y 



w = 



E 



m 



1 + 
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